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A scalar field interacting differently with dark matter and baryons may explain why their ratio is 
of order unity today. We provide three working examples, checking them against the observations 
of CMB, Large Scale Structure, supernovae la, and post-newtonian tests of gravity. Such a scenario 
could make life much easier for supersymmetric dark matter candidates. 
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I. INTRODUCTION 

One of the most debated puzzles in cosmology is why 
the ratio between the energy densities of Dark Energy 
(DE) and Dark Matter (DM) is around two today. Much 
less perceived as a disturbing 'coincidence' is the fact that 
the ratio between the DM and baryonic energy densities 
is around five. 

From an aesthetic point of view, the latter coincidence 
looks more digestible than the former. Indeed, contrary 
to the DE/DM case, the DM to baryon ratio is constant 
in time and its actual value of order unity is not a pecu- 
liarity of the present epoch. Still, there is a priori no rea- 
son why two uncorrelated physical processes, like baryo- 
genesis and the generation of DM abundance, yield relics 
of roughly the same energy densities. 

The purpose of this letter is to present a very sim- 
ple scheme to address this issue: the ratio of the DM to 
baryon energy density is given by a ratio of coupling con- 
stants, like in eq. (21), and so may naturally be of order 
unity. 

In our framework, baryons and DM couple differently 
to the same scalar field, so that the ratio between their 
energy densities, Pb/pd, becomes a field-dependent quan- 
tity. Under general circumstances, the field equations 
of motion exhibit fixed points so that the ratio is driven 
from its early value, determined by baryogenesis and DM 
freeze-out, to a fixed point. 

The idea of a scalar field with different couplings to 
the dark and the visible sectors has a long history, dat- 
ing back to the work of Damour and Gundlach [1], see 
also [2] . More recently, dark matter particles with a time- 
dependent mass, were proposed as a way to increase the 
age of the universe [3], and, after the observation of the 
acceleration of the universe [4], as a way to solve the 
coincidence problem between DE and DM [5]. More in- 
vestigations of new scalar interactions in the dark sector 
were presented in [6,7]. 

Our scenario differs from the previous ones mainly in 
two respects. First, the coupling of the scalar field to 
baryons was generally discarded in the past, both be- 
cause the focus was mainly on the DE/DM sector and 
because it is a potentially dangerous source of fifth forces 
and violations of the equivalence principle. This coupling 



is, on the other hand, crucial here, in order to fix the DM 
to baryon ratio dynamically. Second, in models in which 
DM tracks DE, the growth of cosmological perturbations 
is deeply modified in the past, as a consequence of the 
new scalar interaction and of the exotic dilution rate for 
the DM energy density pdm o.~^ [5]. Here, as we will 
see, the scalar field is fixed in the past and the growth 
of perturbation is just like in the standard ACDM case 
- with at most small deviations in the recent past- re- 
gardless of the strength and non- universality of the scalar 
interaction with DM and baryons. 

Although the main features of our scenario can be im- 
plemented in a generic model, in order to be definite we 
work in the framework of Scalar- Tensor theories of grav- 
ity with purely metric couplings but different conformal 
factors in the visible and the dark sectors, a generaliza- 
tion of the model studied in [1]. 

From a phenomenological point of view. Scalar- Tensor 
theories have the virtue of preserving the equivalence 
principle by construction, i.e. without the need to tune 
the coupling to the visible sector to extremely tiny values. 
From a theoretical point of view they are an attractive 
way to accommodate ultralight scalar fields. Indeed, gen- 
eral covariance implies that the contribution of radiative 
corrections from the (visible and dark) matter sector to 
the scalar field mass is at most of order A'^/M^, A being 
the cosmological constant (see, for instance, [8]). Thus, 
the lightness of the scalar field is just a manifestation 
of the smallness of the cosmological constant, or of the 
curvature of the universe. 

In the following, we will present the general frame- 
work (sect. II), and then give three possible realizations, 
depending on the form of the scalar effective potential 
(sects. Ill, IV and V). Then, in sect. VI, we will make 
some comments on the possible phenomenological impact 
of this scenario. 



II. VISIBLE AND DARK SECTOR 

We will consider an action of the form 

S^Sg+Sb + Sd, (1) 
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where Sg is the gravitational part, given by the sum of 
the Einstein-Hilbert action and the scalar field one, 



(2) 



whereas Sj, and Sg represent the visible (or 'baryonic') 
action and the dark one, respectively, 



Sb,d = Sb,d[^b4^Al^^iip)gf,^] 



(3) 



with "ifb {"ifd) indicating a generic field of the visible 
(dark) sector coupled to the metric A'^{(p)gij,i, {A'^{ip)gfj_i,). 
Since, for our purposes, the most relevant component of 
the visible sector is the baryonic one, we will label fields, 
densities, etc., from this sector with the letter "b". V{(p) 
can be either a true potential or a (Einstein frame) cos- 
mological constant, V{(p) — Vq- 

The emergence of such a structure, with different con- 
formal factors for the visible and the dark sectors can be 
motivated in extra-dimensional models, assuming that 
the two sectors live in different portions of the extra- 
dimensional space. One could also generalize the above 
structure to include more dark matter species, each with 
a different conformal factor, as was done in [7]. Moreover, 
it should be noticed that our results hold for a more gen- 
eral lagrangian than the one above. In practice, as we 
will see, all we need is that the ratio of baryon masses 
to dark matter ones has a non-trivial field dependence. 
In order to preserve the equivalence principle, we have 
to impose that the ratio of proton to neutron masses is 
constant, which is the case for the action in eq. (3). 

If Ab is constant, then the Sg + Sb sector is just that of 
GR, plus a scalar field decoupled from the physical sector. 
The distance to GR is then measured by the quantity 



ab 



dlog Ab 

dip 



(4) 



We consider an homogeneous cosmological space-time 

ds'^ = dt^ -a^{t) df , 

where the matter energy-momentum tensors, T^''^ = 

2{—g)~^^'^5Sb,d/5g'^" admit the perfect-fluid representa- 
tion 

= {Pb,d +Pb,d) '^n'^v - Pb,d g^v , 

with gi^jy u^u" = 1. 

The Fricdmann-Robertson- Walker (FRW) equations 
then take the form 



a 1 



a 6MS 



[pb + 3pb + pd + 3pd + 2M,V^ - 2F] 



(5) 
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Pb + Pd + -^y + ^ 



(6) 



dV 

ctb{pb - 3p5) + ctdiPd - 3pd) + 



(7) 



where is the analogous of eq. (4), and the Bianchi 
identity is 

d{pb,da^) +Pb,d da^ = {pb,d - 3 Pb,d) a^dlog A,d(<p)- 

(8) 

In the following, we will consider only pressureless matter 
both in the visible and in the dark sector, which by eq.(8) 
scale as 



Pb,d ^ Ab,d{y^)a~ 



(9) 



The main point of this paper is that if Ab{ip) ^ Ad{ip) 
the ratio Pb/ Pd becomes a dynamical quantity with a non- 
trivial (^-dependence. Then the present day value of this 
ratio is determined, in principle, both by the early-time 
mechanisms of baryogenesis and freeze-out and by the 
late time behavior of the field (p. 

The above equations have been gven in the so-called 
'Einstein frame'. Actually, an observer using non- 
gravitational 'visible' units, like atomic clocks, wave- 
Icnghts, and masses, is using 'physical (or Jordan) frame' 
units, defined by the rescaled metric 

g^u = Al{ip)g^^ . (10) 

In this frame, pb has the canonical scaling law, pb ^ a 
whereas pd ~ Aci{(f )Ab{(f )~^d~^, the ratio being frame- 
independent, d = aA{ip) is the scale factor in the visible 
frame. In the following, wc will solve the equations of 
motion in the Einstein frame, which is more practical, 
and will then translate the results in the physical one. 

Next, we consider the linear perturbations. We use as 
basis variables the perturbation of the scalar field, Sip, 
and the combination, 



Sd,b = Spd,b/Pd,b - ad,bSip, 



which, since 



Pd,b 
Ad,b 



nd,b, 



(11) 



(12) 



can be identified with the linear perturbation in the 

particle number density. The use of Sd.b instead of 
Spd,b/Pd,b is motivated by the fact that the former are 
(5<y?-independent, as we can read from eq. (9), so that we 
can single out Sip and compute its effective mass properly 
(see eq. (14)). Notice, however, that energy and number 
density perturbations tend to coincide in the Newtonian 
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limit, k^/a^H^ ^ 1 (fc being the comoving momentum), 
in which they differ only by 0{a^ ^a^H^ /k^) terms. In- 
deed, in this limit, the equation for the scalar field per- 
turbation 5^p gives 



fc2 



(13) 



The effective mass m^, has two contributions, 



2 2,2 

m^ = my + nip, 



with, 



1 (fV 



'^V^TT^-TT^ = 1^ f"i + ) • (14) 

In ref. [9] the first term in brackets in to^ was missed 
as a consequence of discussing perturbations in the 
iSpd,b/pd,b, S(p) basis. 

In the Newtonian limit, number density perturbations 
satisfy the the same equations as the energy density per- 
turbations, that is, 

(H' \ 3 

2 + + ad^'j 5'a - -{idd^dSd + Idb^bSb) = , 

S'b + ( 2 + ^ + otbf' j S'b- 2^1bd^d6d + Ibb^bh) = , 

(15) 

where 



7ij = 1 + 2aiaj 1 + 



k^ 



{i,j = b,d), (16) 



from which one can clearly read the extra force due to 
the scalar field, giving a non-universal coupling effective 
at scales such that fc^ > o^m^. This gives a species- 
dependent Yukawa correction to the force between two 
trace pointlike and static masses as 

Fij = ^^{l + 2a,a,e-"^-^) , (17) 

where (i, j — b, d). As we will discuss in sect. IV, if the 

masses arc not pointike, the couplings appearing in (17) 
can be effectively much smaller than the 'bare' values 

The corresponding equations in the physical frame can 
be obtained straightforwardly. Indeed, since pi = A^'^pi 
{i = d, b), we have 



Sb,d = SbM - 4afc5(y9, 



(18) 



so that the two coincide in the Newtonian limit. More- 
over, since rh^a/k = m^a/k, the translation from the 



Einstein equations above to the physical frame simply 
amounts to the change of variable 



d log d = d log a + abdip . 



(19) 



In the following, we will discuss three explicit realiza- 
tions of the dynamical determination of the baryon to 
dark matter ratio, considering three different possibili- 
ties for the scalar potential V{ip). 



III. THE COSMOLOGICAL CONSTANT 

As a first example, we consider the case in which the 
Einstein frame potential is just a constant. Then, if both 
visible and dark matter are pressureless, the right hand 
side of eq. (7) is proportional to the function. 



UbPb + adPd , 



(20) 



where the (/3-dependence resides both in the a's and in 
the p's, sec cq. (9)), and gq = 1. 

The zeroes of the function F{(p) determine the behav- 
ior of the system. If F{ip) has no zeroes, then the field ip 
runs forever and the pb/pd ratio is therefore time depen- 
dent. The Newton constant in the physical frame is also 
time dependent if ab ^ 0, since G/G ~ Habdip/dloga 
and, as discussed in refs. [5,6,10], the non-standard ex- 
pansion rate and non- growth of perturbations may run 
into problem with CMB and large scale structure obser- 
vations. 

We will concentrate on the case in which F(ip) has 
zeroes. Then, it is easy to realize that a zero with a 
positive (negative) slope is an attractive (repulsive) fixed 
point for the field evolution. Depending on the initial 
conditions, the field will proceed towards the next zero 
with a positive slope, which will determine the late-time 
cosmology. In particular, the ratio of densities will be 
given by 



Pb 
Pd 



(Xbifo) ' 



(21) 



ipo being the attractive zero of F{ip) mentioned above. 
The dynamics of the approach to the fixed point depends 
on the form of the function F. In general, the steeper F 
around (po, the faster the field oscillations around it are 
damped. We will discuss this point more explicitly in 
sect. V. On the fixed point the perturbaions exhibit an 
interesting feature. We look for solutions of the equa- 
tions (15) of the form 5d a"^, Sb = bSd- Due to the 
dependence of the "fij on the a's and on fc, the exponent 
m and the bias parameter b are generically expected to be 
ab,d and fc-dependent too. However, due to the relation 
in (21), a 'miracle' happens. We find the solutions 

b=l and m = 1, — 3/2 , 
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c,, -It Jl-48n,,n,,(l + «2„,y^'-)-i 
b = — and m = ; 



(22) 



where we recognize the (fc- independent) m = \, h = 1 
mode typical of a standard matter dominated epoch. 
Thus, independently on the strengths of the scalar in- 
teractions, the stationarity condition (21) implies that 
the fluctuations of the two fluids will always have the 
standard adiabatic growing mode. The other mode cor- 
responds to a potentially dangerous isocurvature pertur- 
bation (Od^d + Q.b5b = 0). However, the growth of this 
mode is extremely suppressed once the post-ncwtonian 
bounds are considered, as we will discuss now. 
The eff'ective mass on the fixed point is 



Mi 



Q.d — ab + 



d\og{ad/ab) 
dip 



(23) 



which is 0{ab^dH'^) ■ Then, unless one of the couplings 
ab_d is enormous, the scalar field is masslcss on solar sys- 
tem scales and the coupling to visible matter is bounded 
by solar system tests of gravity. The recent analysis of 
the delay of radio signals from the Cassini spacecraft [11] 
gives the stringent upper bound 



al < 2.4 X 10"^ (at 3(7). 



(24) 



ad is then fixed by eq. (21), once the best fit value for 
the pb/pd ratio [12], 



Pb_ 
Pd 



= 0.20 ±0.02, 



(25) 



has been imposed. Other astrophysical and cosmological 
bounds on Ud, as discussed for instance by Frieman and 
Gradwhol in [2], would be much less restrictive. 

As anticipated, we conclude that the isocurvature 
mode is likely to play no role, since the bounds above, 
once inserted in eq. (22), imply that the growth exponent 
m is at most O(10~^), to be compared to unity for the 
adiabatic one. 

Thus, once the field has settled to the fixed point, the 
cosmology is just the standard one apart from the value 
of the Pb/ Pd ratio which would be fixed by <^o rather than 
by the standard mechanism for the generation of the relic 
abundance via freeze-out. 



IV. THE CHAMELEON CASE 

Now we consider the case of a quartic potential. 



V = Vq + V- 



4! 



(26) 



Indeed, as we will see shortly, a potential of this type can 
fix the Pb/pd ratio, but cannot account for DE, so that 



a cosmological constant term Vo should be added, 
field attractive fixed point is now given by 



abpb + adPd \ 



1/3 



V 



The 



(27) 



where it is assumed that {abPb + ctdPd)/y ^1- is 
then squeezed very close to the minimum of the poten- 
tial, which fixes the pb/ Pd ratio. The contribution of the 
field-dependent part of the potential to the total energy 
density is suppressed by a large V, 



Pb + Pd 



Pb + Pd 
V 



1/3 



abPb + adPd 
Pb + Pd J 



4/3 



«1, 
(28) 



and, as we anticipated, it cannot account for the observed 
DE if pb^d < V. 

By the same mechanism, the field evolution takes place 
on a timescale much larger than the Hubble one, so that 
we have, for instance 



H 



ab 



dip 



dloga 



f 1 \ 1/3 

tl\otb\\^ ^ 1 <i?, 

(29) 



which can be safely reconciled with existing bounds on 

the variation of G. 

The scalar field gets an effective mass 



m„ 



m 



1 d^V 

Ml dip? 

Pb + Pd f V 



2M2 \pb + pd 



1/3 



abPb + adPd 
Pb + Pd J 



2/3 



(30) 



which, using the value of the cosmological density of mat- 
ter today, gives a range of order 



5 • 10^ Km 



Mi 



-1/6 



(31) 



For ad = 0{1) or larger, the linear perturbations would 

be affected by the scalar field only at scales smaller than 
about 100 Km. For smaller ad the scalar field range 
increases, but at the same time the coupling decreases, 
so that no dramatic effect on the growth of perturbations 
is expected, at least at linear order. 

Since the mass in eq. (31) is much larger than H the ap- 
proach to the minimum of the chameleon potential takes 
place on a time-scale typically much smaller than the age 
of the universe. 

Inside our Galaxy {p ~ 10~^^g/cm^) the range is about 
fifty times smaller than eq. (31). At first sight, such 
macroscopic ranges appear to be in conflict with fifth 
force experiments in the lab or in the solar system. How- 
ever, as it was emphasized in [13] the density-dependence 
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of the effective mass greatly reduces the strength of the 
existing bounds. Indeed, inside baryonic test bodies the 
effective mass is of order 

m-i„ ~ 0.1mm(p6[g/cm3])-i/3 (^^^ " . (32) 

If the size R of the body is much larger than l/m^m 
the field profile is flat inside it except for a shell of size 
~ l/rriipin and the scalar force do not penetrate inside 
it, a phenomenon analogous to the Meissncr effect for 
superconductors, which was dubbed 'thin-shell effect' in 
ref. [13]. In this context, it can be described as an effec- 
tive coupling 

replacing ab in eq. (17) if m^inR ^ 1- Following ref. [13], 
the existing tests on a fifth force translate in an upper 
bound of order 

a,^<O(10), (34) 

relaxing the bound on ai, by many orders of magnitude 
with respect to eq. (24). Geophysical and astrophysical 
searches of deviations from General Relativity give rise 
to much weaker bounds, due to the larger amount of 
suppression caused by the thin-shell effect in the Earth, 
the Moon, and the Sun. 

V. DYNAMICAL DE 

In the previous two sections the role of DE was played 
by a cosmological constant. Here, we examine the pos- 
sibility that the same field ip controlling the pb/pd ratio 
also accounts for a time-dependent DE. In turn, this will 
imply time-dependent Pb/pd ratio, ab, and ad, though 
only for the very last red-shifts. All this might lead to 
conflicts with cosmological observations. 

Indeed, the expansion rate in the physical frame is 
given by, 

H=^Ab{^)-'(l + abJ^). (35) 
a \ a log a / 

We must check that it differs by no more than 10 % from 
the standard one at Big Bang Nucleosynthesis [18], that 
it doesn't shift the CMB peaks by more than 3 %, as 
measured by the 'shift parameter' [14] 

Uo = fll/^Ho r{zdec) = 1-716 ± 0.062, (36) 

fizdec) being the distance to the decoupling surface, and 

that it gives an acceptable age for the universe, -ffo^o = 
0.97 ± 0.01, where we have used the results from WMAP 
at 68% c.l. [15]. 



At small redshifts, the growth of perturbations will de- 
part from the standard solution of a matter-dominated 
epoch, m = b = 1, not only because DE is overtaking 
DM, but also because, when the field moves, the station- 
arity condition of eq. (21) is violated. We will then use 
the constraint from 2dF survey on /3 = f /h = 0.49 ± 0.09 
[17] and b = 1.04±0.11 [19] at the effective redshift 0.15, 
to get 

/(5~ 0.15) = 0.51 ±0.11, (37) 

where / = dlogSd/dloga. 

Moreover, we must check the resulting expansion 
rate at recent epochs against the supcrnovae data [4]. 
We tested the distance-redshift relations of our models 
against the "gold" set of 157 SNe la of ref. [16], using 
flux-averaging statystics as implemented in [20] . 

On top of that, the post-newtonian bound of eq. (24) 
must be satisfied at 5 = 0. 

In order to be be definite, we must chose a functional 
form for the scalar potential and the functions Ab^dif)- 
Of course, lacking a deep theoretical inspiration, a de- 
gree of arbtrariness is present at this stage. However, 
since - after the early-time relaxation to the fixed point 
- the field starts to move again only at small redshifts, its 
late-time behavior mostly depends on the values of the 
functions and their lower derivatives rather than on the 
full functional form. Therefore, the main features of the 
results we will present are independent on the particular 
functional choice. 

Ab{ip) and Ad{ip) must have opposite slopes, in order 
that their logarithmic derivatives satisfy the condition 
(21). Moreover, Ab should be extremely fiat today in or- 
der to meet the post-newtonian bounds of eq. (24). A 
natural way to implement this, is to assume a Z2 sym- 
metry, (f —(p, which sets to zero the first derivatives 
of Ab and Ad at the fixed point = 0. Thus, taking 
into account that without loss of generality we can take 
Ab = Ad = 1 on the fixed point, we can consider the 
expansions 

Ab = l + bbip^ + ■■■ , Ad = l + bd'/ + ■■■ , (38) 

and 

V{^) = Vo + 1/2^^ + Viip'' + ■■■. (39) 

Some illustrative results of solutions passing all the above 
mentioned tests are reported in the figures, together with 
the ACDM case (continuous line). 

In fig. 1 we show the evolution of the scalar field from 
a very early epoch to the present one. As we see, the 
attraction to the fixed point is effective, so that it is 
reached well before nucleosynthesis. Notice that during 
radiation domination the drag towards the fixed point is 
mainly due to the massive particles in the thermal bath. 
Each time one of these becomes non relativistic as the 
temperature decreases, the field receives a 'kick' towards 
the fixed point, as was discussed for instance in ref. [23]. 



5 



Moreover, the trace anomaly gives also a contribution to 
the field evolution, which we have also included in this 
paper. 

As a consequence of the field evolution the pb /pd-ratio 
evolves from its early-time value, set by baryogenesis and 
WIMP decoupling, to the fixed point value of eq. (21). 
Afterwards, it remains practically constant through the 
epochs of nucleosynthesis and matter-radiation decou- 
pling. Finally, depending on the sign of V2, it may or 
may not start a late-time evolution, see fig. 2. 

A positive V2 gives a behavior practically indistinguish- 
able from a pure cosmological constant. On the other 
hand, if V2 is negative, a Z2-breaking phase transition 
takes place when abPb + ctdPd + dV/ dtp becomes negative, 
which happens at small redshifts. The scalar field then 
rolls down from the maximum at </? = and the Pb/ Pd 
ratio changes. The sign of change can be reversed by re- 
versing the signs of hb and hd- The smallness of the effect 
is due mainly to the post-newtonian bound of eq. (24) 
and to the short time available for the field to run. For 
the same reason, the perturbation growth (fig. 4) also 
follows a pattern quite close to ACDM, with / differing 
at most by 10% at redshift less than one. 

The equation of state (fig. 5) of the scalar field could 
be as high as = —0.7 and still pass the SNe test at 95 
% c.L. 
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FIG. 2. The ratio between baryons and dark matter en- 
ergy densities for three illustrative examples satisfying all the 
experimental bounds discussed in the text. The solid line 
corresponds to the ACDM model with Q.a = 0.7. The dot- 
ted ( V2 = -65Vo, 14 = Vo, bb = -0.05, fed = 5), dashed 
( V2 = -lOOVo, 14 = lO^Vo, bb = -0.01, fed = 5) and 
dot-dashed ( V2 = -SOOVb, 14 = 2-lO^Vb, feb = -0.05, fed = 5) 
lines correspond to models passing all the tests mentioned in 
the text. In each case Vo is adjusted to give the right amount 
of DE today. 




Log(z+ 1) 

FIG. 1. The evolution of the field at early time with differ- 
ent initial conditions. The fixed point is reached well before 
nucleosynthesis. In this example we used bb = —0.05, fed = 5. 
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FIG. 3. The coupling functions ab and otd in the examples 
of fig. 2 
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FIG. 4. Evolution of the growth exponent of the linear per- 
turbations, f. 




FIG. 5. Equation of state for the scalar field. 



VI. CONCLUSIONS 

The most popular DM candidate is the lighest super- 
symmetric particle, which, being stable and weakly in- 
teracting, has all the good qualities to be a successful 
WIMP. However, in practice, experimental constraints 
restrict the allowed parameter space of the Minimal Su- 
persymmetric Standard Model (MSSM) to regions where 
the neutralino relic density is typically one order of mag- 
nitude below the WMAP preferred range (see, for in- 
stance, ref. [21]). Bino-like neutralinos, as preferred in 
the mSugra subclass of MSSM, suffer from the oppo- 
site problem, being usually overproduced because of their 
lower annihilation rate. As a result, the allowed param- 
eter space of MSSM after the WMAP results has shrunk 



to regions were extra annihilation channels are present, 
which requires an increasing degree of fine tuning. 

The scenario presented in this paper has the poten- 
tial to alleviate this problem considerably. The DM relic 
abundance produced at neutralino freeze-out, can be 'ad- 
justed' by the late time dynamical mechanism presented 
here, thus rescuing wide portions of the parameter space. 
Of course, once the coupling functions are chosen, the 
fixed point attraction should be strong enough as to drive 
the initial DM abundance to the WMAP value before 
nucleosynthesis, which would provide the limits of the 
region of the rescuable parameter space. 

Other post-freeze-out adjustment mechanisms were 
proposed in [22] and in [23]. In both cases the adjust- 
ment could only be an enhancement of the relic density, 
whereas here we could have a reduction as well, a case of 
interest for a bino-like neutralino. 

An obvious signature of these adjustment mechanism 
would be a laboratory discovery of a DM particle whose 
interactions do not match the relic abundance computed 
in the standard way. 

On the cosmological/gravitational side, an improve- 
ment of the post-newtonian bounds and an increase in 
the sensitivity to the linear growth rate / by at least a 
factor of two might constrain this scenario considerably. 
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